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Abstract 

Statistical mechanics can only be ultimately justified in terms of microscopic dynam- 
ics (classical, quantum, relativistic, or any other). It is known that Boltzmann-Gibbs 
statistics is based on the hypothesis of exponential sensitivity to the initial condi- 
tions, mixing and ergodicity in Gibbs F-space. What are the corresponding hypoth- 
esis for nonextensive statistical mechanics? A scenario for answering such question 
is advanced, which naturally includes the a priori determination of the entropic 
index q, as well as its cause and manifestations, for say many-body Hamiltonian 
systems, in (i) sensitivity to the initial conditions in Gibbs F-space, (ii) relaxation 
of macroscopic quantities towards their values in anomalous stationary states that 
differ from the usual thermal equilibrium (e.g., in some classes of metastable or 
quasi-stationary states), and (iii) energy distribution in the T-space for the same 
anomalous stationary states. 
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1 INTRODUCTION 



There is no reasonable doubt that it was, at least intuitively, neat and clear 
for Ludwig Boltzmann and for Albert Einstein [1,2] that statistical mechanics 
descends, in one way or another, from microscopic dynamics. It should be so 
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for all essential concepts of the theory, including the celebrated expression for 
the entropy, namely (written here for a discrete case) 



w 



SbG ^ -kJ^Pi^^Pi ^ (1) 



from now on referred to as the Boltzmann-Gibbs (BG) entropy. However, and 
interestingly enough, it is until now unknown how to start from say Newton's 
law F = ma and end up with this remarkable functional form and its as- 
sociated variational principle. It is Boltzmann's magnificent intuition which 
made him to propose this specific microscopic expression for Clausius' ther- 
modynamic entropy. Expression which is totally consistent with Newton's law, 
although rigorously speaking we still do not know why. 

If such is the situation for Sbg, it is clear that the situation has to be even 
more intriguing - to say the least! - for any generalization of Sbg, such as [3] 

S^^k {qeR; Si^Sbg). (2) 



Even if we do not yet know how Sbg descends from dynamics and what re- 
strictive premises the microscopic dynamics must satisfy, it has been profusely 
verified that Sbg is the natural entropic form everytime we have nonintegrable 
and sufficiently chaotic dynamics (i.e., positive Lyapunov exponents), known 
to yield quick mixing and eventually ergodicity in phase space. A natural 
question arises: what happens if the system is nonintegrable but all Lyapunov 
exponents vanish? More precisely, what happens when the sensitivity to the 
initial conditions diverges less than exponentially? It seems plausible to imag- 
ine that the answer depends on how the sensitivity diverges asymptotically, a 
power-law?, a slow logarithmic divergence?, some other type? The case that 
we focus on in the present paper, and generically in nonextensive statisti- 
cal mechanics, is the ubiquitous power-law. It is for such dynamics that one 
imagines that Sq would be the appropriate physical entropy to be used for var- 
ious purposes, including the bridging to thermodynamics. Furthermore, one 
imagines that the appropriate value of the entropic index q should reflect the 
exponent of the power-law, thus constituting the basis for universality classes 
of nonextensivity. 



2 A POSSIBLE TRIANGLE FOR THE ENTROPIC INDEX q 



Although most of the notions that we shall use here occur in all types of 
dynamics, including quantum ones, we shall focus on classical systems for 
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simplicity. 



2.1 SENSITIVITY TO THE INITIAL CONDITIONS 



Let us start witii tiie notion of sensitivity to the initial conditions. At the 
present point it is enough, as an illustration, to consider a one-dimensional 
nonlinear dynamical system (e.g., an unimodal map like the logistic one) 
characterized by x{t). We define the sensitivity (to the initial conditions) 
^ = limAa;(o)^o ^^(^)/^3;(0), where Ax{t) = x{t)—x'{t) denotes the difference 
between two trajectories x(t) and x'{t) that are initially close to each other. 
The following differential equation is typically satisfied: 

§ - • 

where Ai is the Lyapunov exponent (or essentially the largest one in more 
general situations) . The solution of course is 

eW = e^^*, (4) 



hence In ^ oc t . This solution generically holds everytime Ai ^ 0. But if we 
have Ai = (referred to as weak chaos if there is mixing, and integrability 
otherwise), we expect the adequate differential equation to typically be the 
following generahzation of Eq. (3): 



dt 



K.ni'- , (5) 



where sen stands for sensitivity. Its solution is given by 

eW = eJ-«*, (6) 

hence ln,^^„ ^ oc i , with [4] 

= [1 + (1 - g)a;]V(i-«) (e^ = e") , (7) 



and 

X^-Q — 1 

luq X = — (Ini X — Inx) . (8) 

(We remind that is always real and nonnegative; for q < 1, it vanishes for 
all values of x such that 1 + (1 — q)x < 0). 
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If Ai = and, in spite of that, there is mixing in the full dynamical space 
of the system, we typically expect g^g^ < 1 and Ag^^„ > 0; if there is no 
mixing, we typically expect qsen > 1 and \q^^^ < 0. It is of course possible in 
even more pathological cases (and therefore more rare in natural or artificial 
systems) that both \\ and Ag,^„ vanish, but such cases are out of the scope 
of the present discussion, centered on power-laws (and not on even weaker 
behaviors such as the logarithmic ones). 

It is by now well known that an illustration of the situation we are focusing 
on here is provided by the edge of chaos of the logistic map: qsen = 0.2445... 
and \q^^„ = 1/(1 — qsen) = 1.3236... [5,6]. Many other paradigmatic dynamical 
systems arc known or expected to also enter, in one way or another, into the 
same category. This is the case of globally coupled standard maps close to 
integr ability ([7,8,9] and references therein) or long-range-interacting many- 
body Hamiltonian systems ([10] and references therein). 

Let us finally mention that, at least for simple one-dimensional maps, qsen 
also enters prominently in basic properties such as the entropy production per 
unit time [11] (essentially the Kolmogorov-Sinai entropy rate), and the values 
at which the multifractal function (associated with the edge of chaos) may 
vanish [12]. 



2.2 RELAXATION 



We define the quantity 



where O is a macroscopic observable relaxing towards its value at a stationary 
state (thermal equilibrium, or some other). The associated differential equa- 
tion very frequently is 



Ti > being the relaxation time. The solution of this equation is 

Q{t) = e-*/^i . (11) 



The quantity ri is in principle expected to depend on the observable O. How- 
ever, in all generic cases one expects (according to Krylov's classical remarks) 
1/ti ~ Ai if Ai > 0. This is quite natural since the relaxation of a macroscopic 
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quantity reflects nothing but the exploration that the system makes of its fuU 
dynamical space in its search of a possible stationary state (given some class 
of initial conditions), and this occurs at a rythm dictated by the Lyapunov ex- 
ponent Ai (or, in general, by the full Lyapunov spectrum). In particular, if Ai 
(or the maximal Lyapunov exponent) vanishes, we expect ti to diverge, thus 
making Eq. (11) inappropriate. It is quite plausible that, in such situation, 
Eq. (10) becomes generalized into 

^ = — -n^^^^ (12) 

at Tg , 



where rel stands for relaxation. The solution is of course given by 

m = e',:i:'-^ (r,.,>0) (13) 



being typically g^e/ > 1- It could well happen that g^e/ > 1 depends on the 
specific observable O we might be interested in, but even if such is the case, we 
expect all possible qreis to be simply inter-related. Moreover, we expect them 
to be essentially characterized by the speed at which the nonlinear dynamics 
makes the system to approach its stable (or metastable) stationary (or quasis- 
tationary) state, given the class of initial conditions within which the system 
has been started. More specifically, when Ai = 0, we expect simple hypothesis 
such as ergodicity in full phase space (basic for the validity of BG statistical 
mechanics) to be not fulflUed. Even more, we expect (under conditions to be 
determined) the system to have a tendency of lengthily living in some kind 
of attractor or pseudo-attractor, whose Lebesgue measure in full dynamical 
space would typically be zero, as it is the case of (multi) fractals or scale-free 
networks. If we assume an ensemble whose Lebesgue measure is different from 
zero ai t — inside a special region of full space, we might observe a shrink- 
ing of the Lebesgue measure as time goes on. The loss of Lebesgue measure 
would possibly occur through a ^-exponential function whose value of q would 
be noted qj-ei- From now on, we reserve the notation qrei for the shrinking of 
Lebesgue measure, if such is the case. All other g^ez's would be simply related 
to this one. Let us emphasize that the present scenario does not exclude that, 
at even larger times, the system escapes from this (multi) fractal "prison" and 
ultimately satifies ergodicity, thus making a crossover from g^e/ to g = 1. This 
would in fact be the case unless some limiting situations (see [8]) are taken 
which strictly guarantee that Ai vanishes. When such escape is possible, aging 
phenomena might be present [13,14], aging being presumably related to a slow 
collective dynamics which would prepare the escape from the "prison" . 

A typical illustration of the picture we have described in the present subsection 
is, as before, the edge of chaos of the logistic map, where it is qrei — 2.4 [15]. 
Other paradigmatic illustrations should include the systems mentioned in the 
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previous subsection, but this remains to be verified. 

Let us finally mention that, at least for simple one- dimensional maps, qrei also 
emerges [15] in basic properties such as the fractal dimension (maximum of 
the multifractal function) of the attractor at the edge of chaos. 



2.3 STATIONARY STATE 



Let us now address what may in some sense be considered as the most useful 
thermostatistical function, namely the distribution of total energies in Gibbs 
F-space of a thermodynamical system in contact with a canonical thermostat 
fixing the temperature T. 

Let us first consider a Hamiltonian system whose elements interact only locally, 
or almost not at all (e.g., the ideal gas). A typical example of local interaction 
is a classical two-body (attractive) potential which has no singularity at the 
origin, and which asymptotically decays with distance r as — 1/r" with a/d > 
1, where d is the space dimension of the system; for instance, for the usual 
Lennard- Jones fluid, we have a = 6 and d = 3. We know that the relevant 
stationary state is thermal equilibrium and that 



where {Ei} is the set of eigenvalues of the Hamiltonian with specific boundary 
conditions. This celebrated distribution can be obtained through the textbook 
procedure which consists in optimizing Sbg under the constraints J2iPi = 1 
and Y^iPiEi — Ui, where Ui is the total internal energy. It seems to be realized 
by only a few physicists that this distribution can also be obtained, albeit 
without formal justification up to now, through a differential equation, namely 

^ = -A(P.^.). (15) 



The black-body radiation law was first found, and published in October 1900 
[16], by Max Planck through an heuristic differential equation, precisely one 
which, with appropriate variables, contains Eq. (15) as a particular case! This 
is known only by those interested in the history of science. And it has not 
been incorporated in standard textbooks presumably because of the lack of 
formal justification. Nevertheless the mathematical fact is there, and in my 
opinion it is perhaps not a mere coincidence. Indeed, since the sensitivity to the 
initial conditions follows an exponential law, and the same typically does its 
prominent consequence namely the relaxation of macroscopic thermostatistical 
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quantities, it cannot be a very great STirprise that the same does happen to 
the distribution associated with the stationary state towards which the system 
relaxes. 



For Hamiltonian systems whose elements do not interact locally but globally 
(e.g., the classical systems mentioned previously but with < a/d < 1, in 
which case the potential is not integrable), it seems quite plausible that the 
Eq. (15) becomes generalized into 

— PqstaAPl'^qstat) ) l-LOJ 



where stat stands for stationary state, and /^^^^^^ is the adequate inverse tem- 
perature (in units of k). The solution is given by 

Pi = — {^..t ^ ykT^taf, ^ E ^iSf , (17) 

Qstat j 



and one typically expects qstat > 1. It is in fact precisely this solution which 
emerges [3] by optimizing, through conveniently written norm and energy 
constraints, Sg instead of Sbg- 

There is molecular dynamical numerical evidence [17] for some classical long- 
range-interacting many-body inertial rotators on (i-dimensional lattices that 
they might satisfy, in some robust metastable states in the presence of a ther- 
mostat, the total energy distribution indicated in Eq. (17). But, even at the 
numerical level, the problem is uneasy, though not impossible, to study. This 
is the subject of ongoing efforts. The strategy which is followed is to consider 
an isolated system of N rotators at a fixed total energy. Then consider once for 
ever a subset of M rotators among the available. Finally, one must numeri- 
cally approach the t oo limit of the M ^ oo limit of the N ^ oo limit. The 
whole procedure must be repeated many times {N realizations ~^ mathemat- 
ically speaking) for realizations differing in the precise multi-particle initial 
conditions, always within some specific and nontrivial class of initial condi- 
tions, and then averaging the results. The distribution obtained in this manner 
(i.e., limt^oo liniM^oo liniAr^oo liniAr^^^;.^^^.^^^-^oo) is to be compared with Eq. 
(17). 

Let us make one more remark, which concerns the distributions of momenta 
of these rotators. If we denote with EI({Lj}, {9i}) the Hamiltonian associated 
with say planar rotators, the distribution in F-space is given by 

p({L,},{^,})«e-f--««^'M^^». (18) 
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Then the one-body angular momentum marginal distribution p{L) is given by 



We do not yet know what is the result for this distribution but it could also 
be a g'L-exponential in with ql > 1 related to Qstat is some still unknown 
manner. One expects (essentially because the integral of a g-exponential is 
another g-exponential) the index to either coincide with qstat or be a simple 
function of it (sec [18]). In any case, the available numerical evidence is not in- 
consistent with such posibilities ([8] and references therein). See, for instance. 
Fig. 1 (from [8]), where for the first time the momentum distribution is calcu- 
lated (through molecular dynamics) for a system which is sensibly closer to a 
canonical ensemble (contact with a large thermostat) than the microcanonical 
ensemble considered up to now in the literature for such systems. 



The g-indices discussed in Section 2 and their possible interconnections are 
schematically represented in Fig. 2. It is our scenario that they are differ- 
ent "faces" of the same basic phenomenon, namely how the system "hkes" to 
evolve and mix in part or all of its full phase space (Gibbs F-space for Hamil- 
tonian systems) given its initial conditions. If the system is strongly chaotic 
(i.e., positive Lyapunov exponents), essentially satisfying Boltzmann's "molec- 
ular chaos hypothesis'' , then the system will be mixing and ergodic allover 
the entire phase space. BG statistical mechanics is expected to appropriately 
describe the thermal statistics of the system, the associated microscopic ex- 
pression for the entropy undoubtedly being Sbg- But if the system is only 
weakly chaotic (i.e., zero Lyapunov exponents, but not integrable), then the 
occupancy of phase space might be much more complex everytime the system 
is started within specific nonzero- Lebesgue-measure regions of phase space. 
The region of phase space where the system may lengthily evolve before even- 
tually occupying the entire space might have a (multi)fractal or hierarchical 
structure of the type currently referred to as scale-free network (see [19] and 
references therein). As long as the system remains in that region, its ther- 
mostatistics might be adequately be described within nonextensive statistical 
mechanics, the associated microscopic expression for the thermodynamic en- 
tropy presumably being Sq with a specific value of q. During the metastable 
or quasi- stationary state we would then have {qsen,(lrei, Qstat) 7^ (1, 1, 1), and 
at later times (in many cases, infinitely later!) a crossover would occur to 
Qstat — Qrei — Qsen — 1- This Confluence onto a single value of q for BG 



piL) = I (Ul, dU) {Ul, de,) p({L,}, {6,}) 




(19) 



3 FINAL REMARKS 
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Fig. 1. Snapshots of the evolution of PDF for momenta in the (L, 0) — space (i.e., in 
^-space), starting with a "liot water bag". In solid black circles, the instantaneous 
PDF at time t = lO'' (with A; = 0, 1, ...7), for M = 500 spins inside an = 5000 
spin system. In empty circles, the average of 10^ realizations for TV = 10^ spins (all 
averages made during the QSS plateau). In dashed line, the g^-exponential fitting 
curve with T = Tqss = 0.38 and ql = 3.7. Finally, in solid line, the analytical 
Gaussian PDF. These last three curves are the same in every frame, and are plotted 
for reference. Prom [8]. 
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Fig. 2. The triangle of the basic values of q, namely those associated with sensitiv- 
ity to the initial conditions, relaxation and stationary state. For the most relevant 
situations we expect qsen < l, Irei ^ 1 and qstat > 1- These indices are presumably 
inter-related since they all descend from the particular dynamical exploration that 
the system does of its full phase space. For example, for long-range Hamiltonian 
systems characterized by the decay exponent a and the dimension d, it could be 
that qstat decreases from a value above unity (e.g., 2 or 3/2) to unity when a/d 
increases from zero to unity. For such systems one expects relations like the (par- 
ticularly simple) qstat = Irei = 2 — qsen or similar ones. In any case, it is clear 
that, for a/d > 1 (i.e., when BG statistics is known to be the correct one), one 
has qstat = Qrei = Qsen = 1- All the weakly chaotic systems focused on here are 
expected to have well defined values for qsen and qrei, but only those associated 
with a Hamiltonian are expected to also have a well defined value for qstat- 

statistics might be at the origin of occasional confusion within the interested 
community. Indeed, three essentially different (though related) concepts just 
happen to coincide within the BG formalism. 
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